^ 



TUW-04-14 

Collective modes of an Anisotropic Quark-Gluon Plasma II 

Paul Romatschke and Michael Strickland 

Institut fiir Theoretische Physik, Technische Universitdt Wien, 
Wiedner Hauptstrasse 8-10, A-IO4O Vienna, Austria 



Abstract 

We continue our exploration of the collective modes of an anisotropic quark gluon plasma by 
CD ' extending our previous analysis to arbitrary Riemann sheets. We demonstrate that in the presence 

^^^ . of momentum-space anisotropics in the parton distribution functions there are new relevant singu- 

^ , larities on the neighboring unphysical sheets. We then show that for sufficiently strong anisotropics 

^ I that these singularities move into the region of spacelike momentum and their effect can extend 

down to the physical sheet. In order to demonstrate this explicitly we consider the polarization 
tensor for gluons propagating parallel to the anisotropy direction. We derive analytic expressions 
for the gluon structure functions in this case and then analytically continue them to unphysical 
,.^ , Riemann sheets. Using the resulting analytic continuations we numerically determine the position 

QQ I of the unphysical singularities. We then show that in the limit of infinite contraction of the dis- 

Q^ ' tribution function along the anisotropy direction that the unphysical singularities move onto the 

\^ • physical sheet and result in real spacelike modes at large momenta for all "out-of-plane" angles of 

O ■ propagation. 

o 



o 






X 



PACS numbers: 11.15Bt, 04.25.Nx, ll.lOWx, 12.38Mh 



I. INTRODUCTION 

The ultrarelativistic heavy ion colHsion experiments ongoing at the Brookhaven Relativistic 
Heavy Ion Colhder (RHIC) and planned at the CERN Large Hadron Colhder (LHC) will study the 
behavior of nuclear matter under extreme conditions. Specifically, these experiments will explore 
the QCD phase diagram at large temperatures and small quark chemical potentials. Based on 
the data currently available from the RHIC collisions it seems that a thermalized state has been 
created during the collisions [l|. Remarkably it seems that the thermalization proceeds rather 
rapidly in contradiction to estimates from leading order equilibrium perturbation theory. However, 
to truly understand how the plasma evolves and thermalizes one has to go beyond the equilibrium 
description. In this paper we expand upon our previous studies of the collective modes of a quark 
gluon plasma which is (at least approximately) homogeneous and stationary but anisotropic in 
momentum space. 

These types of distribution functions are relevant because of the approximate longitudinal boost 
invariance in the central rapidity region of ultrarelativistic heavy ion collisions. This implies that 
the initial distribution functions for the partons are practically delta functions in longitudinal 
momentum. Such an anisotropic quark-gluon plasma appears to be qualitatively different from 
the isotropic one since the quasiparticle collective modes can then be unstable [j, y, Ij, ly, |y, lDj 
la y, [10]. The presence of these instabilities can dramatically influence the system's evolution 
leading, in particular, to its faster equilibration and isotropization. Treating this problem in all 
of its generality is a daunting task. In order to make progress we consider the limit of very high 
transverse temperatures at which the non-equilibrium collective behavior is describable in terms of 
processes in which all loop momenta are hard. In the case of thermal equilibrium hard corresponds 
to momenta of order T but in the non-equilibrium case hard corresponds to an unspecified scale 
contained in the distribution function. 

In a previous paper [3| we calculated the hard-loop gluon polarization tensor in the case that 
the momentum space anisotropy is obtained from an isotropic distribution by the rescaling of one 
direction in momentum space. The resulting expression for the gluon polarization tensor was then 
decomposed into a four-component tensor basis and the structure functions associated with this 
tensor basis where computed numerically for general anisotropics and analytically in the limit 
of small anisotropics. We demonstrated that a contraction of an isotropic distribution function 
along the anisotropy direction, n, resulted in one additional stable quasiparticle mode, two damped 
quasiparticle modes in the lower half plane, and two unstable (anti-damped) quasiparticle modes in 
the upper half plane. In the case that the isotropic distribution was stretched along the anisotropy 
direction we found that again an additional stable quasiparticle mode was generated but only one 
damped and one anti-damped quasiparticle mode were found in this case. The unstable modes 
found correspond to electric or ma gnet ic t ype instabilities with the latter being analogous to the 
Weibel instability in QED plasmas [ij, 0,0,113,0 • 



In this paper we continue our exploration of the collective modes of an anisotropic quark gluon 
plasma by extending our previous analysis to arbitrary Riemann sheets. We demonstrate that in 
the presence of momentum-space anisotropics in the parton distribution functions there are new 
relevant singularities on the neighboring unphysical sheets. We then show that for sufficiently 
strong anisotropics that these singularities move into the region of spacelike momentum and their 
effect can extend down to the physical sheet. In order to demonstrate this explicitly we consider the 
polarization tensor for gluons propagating parallel to the anisotropy direction. We derive analytic 
expressions for the gluon structure functions in this case and then analytically continue them to 



unphysical Riemann sheets. Using the resulting analytic continuations we numerically determine 
the position of the unphysical singularities. We then show that in the limit of infinite contraction 
of the distribution function along the anisotropy direction that the unphysical singularities move 
onto the physical sheet and result in real spacelike modes at large momenta for all "out-of-plane" 
angles of propagation. 

The organization of the paper is as follows: In Sec. ^ we first review the necessary integral 
expressions for the hard-loop gluon polarization tensor and then in Sec. IIIII we present analytic 
expressions for the hard-loop gluon polarization tensor structure functions in the case that the 
gluon is propagating parallel to the anisotropy direction. In Sees. 1111 AinTTOl we extend these 
expressions to arbitrary Riemann sheets and solve the dispersion relations for the singularities 
existing on the neighboring unphysical sheets. In Sec. IIVI we present analytic expressions for the 
gluon polarization tensor structure functions in the large-anisotropy limit and solve the resulting 
dispersion relations for arbitrary angle of propagation. In Sec. we summarize the results and 
speculate about the possible impact of the now relevant unphysical singularities. 

II. GLUON POLARIZATION TENSOR REVISITED 

The hard-loop gluon polarization tensor of an anisotropic system is given by 

where K = (cjjk), V = (l,v), v = p/|p| and 

/(p) ^ 2Nf (n(p) + n(p)) + 4iV,n,(p) . (2) 

Note that in (^ we have specialized to spacelike Lorentz indices; however, it is possible to derive 
the polarization tensor also for arbitrary Lorentz indices. 

To simplify the calculation we follow Ref. [3] and require the distribution function /(p) to be 
given by 



Here /iso is an arbitrary isotropic distribution function, n is the direction of the anisotropy, ^ > — 1 
is a parameter reflecting the strength of the anisotropy, and A^(^) is a normalization constant. 
To fix N{^) we require that the number density to be the same both for isotropic and arbitrary 
anisotropic systems, 

/^ /iso(p) = J^ Mp) = NiO J^ ho (v/p^ + C(p-n)2) ^ (4) 

and can be evaluated to be 

N{0 = VlTI . (5) 

Using an appropriate tensor basis J2| one can then decompose the self-energy into four structure 
functions a, /?, 7, and 5 by taking the contractions 

k'U'^k^ = k'^P , 
n'li'^y = h^k^6 , 

Trrf^' = 2a + /? + 7, (6) 



where n* = (J*-' — k^y /k'^)h^ . The structure functions then depend on u, k and the angle k • n = 
cos 9n as well as on the strength of the anisotropy, ^. Integral expressions for a,f3,^, and 6 for 
arbitrary angle of propagation and anisotropy parameter ^ can be found in Ref. J3]-^ 

III. SPECIAL CASE I : k || n 

Let us now consider the case where the momentum k is in the direction of the anisotropy, fi, 
i.e. 9n = 0. Using the changes of variables 



p'=p^{i+a^-hf) , 



allows us to simplify Eq. (^ to 



where 
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(8) 
(9) 



Taking the contractions in Eq. (jB)) the structure functions can further be simplified using 

k-v = A;n-v = /c cos 9 . (10) 

While one integration becomes straightforward, the remaining integration can be performed after 
a little bit of algebra and one obtains for the relevant contractions 
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z + l + ie 
z — 1 + ie 



, (11) 



where z = oj/k and 6 = 6k. Note that for this angle of propagation the structure function 7 
vanishes. 



A. Extension to unphysical sheets 

The above structure functions all possess a logarithmic cut running along the real z axis for 
z^ < 1. It is, however, possible to extend their definition beyond this cut, which then corresponds 



^ In Ref. [7| N{^) was fixed to be -/V(^) = 1 so the reader should make sure to adjust for the difference 
where appropriate. 



to an unphysical Riemann sheet of z. More precisely, there are two such unphysical sheets which 
can be accessed by either extending the physical sheet from above or below the cut. In order 
to make this extension we first compute the structure functions on the physical sheet using the 
integral representation given in Eq. (jSJ. We can then deform the integration contour which runs 
from cos = —\ along the real axis to cos d = \ into the lower half plane (LHP) and then move the 
point we are interested in from the upper to the lower half plane. Deforming the original contour 
so that it again runs along the real axis we see that we pick up an extra contribution corresponding 
to the residue at the point in the LHP. 

This procedure can be used to analytically continue the structure functions for all values of On 
but in the case that ^„ = the residue can be evaluated straightforwardly and is simply —2Tri. 
The structure functions can likewise be extended into the upper half plane (UHP) and the residue 
is then +27rz. The resulting rule is then the expected one, namely, that the structure functions can 
be extended to unphysical sheets by the usual extension of the logarithm 



In ( ^— y j = In f ^—j j + i arg l^—^ ) + 27rn 



(12) 



where n specifies the sheet number. Note that n = 1 extends the physical logarithm into the UHP 
and n = — 1 extends it into the LHP. Higher n correspond to higher sheets which can be safely 
ignored as we will discuss below. 

B. Collective modes 

The dispersion relations for the gluonic modes in an anisotropic quark-gluon plasma are in 
general given by the zeros of 

A;^^ = k^ -Lo^ + a = , 
A^i = (A:^ - ^2 ^ Q + 7)(/3 - a;2) - k^h'^d^ = . (13) 

In the case k| |n, however, 7 vanishes identically, as does n^ = 1 — (k • n)^. Therefore, it is sufficient 
to solve the equations 

k'^ -uj"^ + a = , 

P-u;^ = 0, (14) 

which will be referred to as a- and /3-modes, respectively. These modes correspond to poles in the 
propagator for a- and /3- modes. We will use this term to describe solutions on both physical and 
unphysical Riemann sheets; however, whenever we are speaking about the unphysical singularities 
we will always explicitly label them as unphysical a- and /3-modes. The reader should be aware 
that these unphysical singularities do not correspond to real degrees of freedom unless the solution 
associated with them moves onto the physical Riemann sheet. 

Before we present the dispersion relations we would like to first count the number of modes on 
the various sheets in order to be assured that we have indeed found all solutions. The number of 
modes can be counted by doing a so-called Nyquist analysis, based on the special case of Cauchy's 
integral, 

where N and P are the number of zeros and poles of f{z) times their multiplicity in the region 
encircled by the closed path C. 
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FIG. 1: Contour C in the complex z plane 
used for finite ^ Nyquist analysis. 



FIG. 2: Contour Clhp in the complex z 
plane used for finite ^ Nyquist analysis on 
the LHP unphysical sheet. 



Choosing f[z) = fe^(l 



+ a{z) and using the explicit form of the structure function given 



in Eq. Hll|l. one finds that f{z) has a logarithmic cut for real z^ < 1, while being analytic for all 
other finite z. One can then choose the contour C depicted in Fig. ^ so that P = 0. We then 
evaluate the respective pieces of the contour C, finding 



2m 



47ri + + + 47riG( 



lini/(z)) 

2^0 



2 + 2e(- 



lim/(z)). 

Z-+0 



(16) 



The first contribution comes from the large circle at \z\ ^ 1, while the first zero is the contribution 
from the path connecting the large circle with the contour around z"^ < 1; the second zero is the 
contribution from the small half-circles around z = ±1. The last contribution comes from the 
straight lines running infinitesimally above and below the cut at z^ < 1. They can be evaluated 
by using 



/'(^) 



l+ie f{z) 



f{l + ie) 



and therefore represent the "winding number" of the image of f{z) around the origin. Since from 
Eq. (|11() it is clear that Re /(—I + ie) > and for z^ < 1 one has Im/(z) = only for 2 = 0, the 
winding number can either be zero or one, depending on the sign of Wvn.z^Q'R.e f {z + ie). For the 
path with Imz < one proceeds similarly, finding the result given in Eq. H16|). 
For the /5-mode, the same techniques can be applied to find the result 



N, 



f3,phys 



(18) 



These two modes, together with two modes coming from Na^phys: correspond to the propagating 
and stable modes for positive and negative (real) frequencies, while the modes that depend on 
the sign of the static limit of f{z) correspond to solutions withpurely imaginary z which are the 
damed and anti-damped physical modes already found in Ref. J7|. 

We now want to extend the above analysis to the unphysical sheets that can be accessed by 
extending the structure functions through the cut to the LHP and UHP, as discussed above. Let 
us first discuss the a-mode in the unphysical LHP by introducing fLHp{z) = k'^{\ — z^) + aLHp{z) 
and choosing a contour Clhp which is shown in Fig. [21 One notable difference to the physical 
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FIG. 3: Dispersion relations for the unphysical a-modes for 0„ = and (a) ^ = 0.1, (b) ^ = 1, 
and (c) ^ = 10. Dotted lines correspond to the modes which depend on the sign of /lhp (type-1). 
Note that in (a), (b), and (c) the type-1 unphysical o-mode only exists for k > 0.184, k > 0.595, 
and k > 1.95, respectively. Dashed lines correspond to modes which do not depend on the sign of 
/lhp (type-2). Solid line in (a) is the isotropic result. Note that the scale changes in each plot. 



sheet is that now fLHp{z) has a pole of second order at z = —i/y/C, so that Pa,LHP 
than that the analysis is conducted as in the case above, finding 



N, 



a,LHP 



°a,LHP = -6(- lim fLHp{z)) , 



2. Other 



(19) 



so that Na,LHP = 1 + 9(lim2_»o fLHp{z)). 

For the UHP, one uses a similar analysis to find 



N, 



a,UHP 



l + e(lim/f/Hp(^)), 



(20) 



while for the unphysical /3-mode one has Np^iHP = ^p^UHP = 2. Closer inspection of these modes 
shows that the unphysical a-modes correspond to solutions with purely imaginary z, while the 
unphysical /?- modes are solutions with complex z with positive/negative real and imaginary parts. 

A plot of the position of the unphysical a-modes for ^ = {0.1, 1, 10} is shown in Fig. |31 Since 
there are two unphysical a-modes we will refer to the ones that depend on the sign of Jlhp as 
type-1 and those which don't as type-2. We can observe from this Figure that, as ^ goes to zero, 
the type-1 mode (dotted line) becomes degenerate with the isotropic solution (solid line in Fig. EJi) 
and the type-2 mode (dashed line) moves off to infinity. In the opposite limit, .^ ^ 1, we find that 
the type-2 mode moves towards the real axis while the threshold for existence of the type-1 mode 
moves to infinity. 

Note that even in the isotropic limit there are unphysical a-modes on the neighboring Riemann 
sheets: one in the UHP and one in the LHP. These unphysical modes are related to the presence of 
dynamical screening of the magnetic interaction. In fact, in the limit of small momentum it can be 
shown that these solutions are directly related to the small z behavior of the isotropic transverse 
polarization tensor. This can be seen by taking the small z and E, limits of ajjHP 



lim lim 



auHP 



f(i.|«-l..o(e 



2 Z^) 



(21) 

The first term in (|21|) comes directly from the logarithmic cut. Taking the isotropic limit, ^ = 0, 
inserting this into the a-mode equation given in Eq. H14|). and then taking the static limit we see 
that there is a solution at purely imaginary uj 

u; = -Xr + OiC,z'). (22) 
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FIG. 4: Dispersion relations for 0„ = 0, ^ = 
10: besides two propagating a- and /3-modes 
on the physical sheet (dashed and full lines, 
respectively) there is also an unphysical /?- 
mode (dotted line) for which we plot only 
the real part of the pole position. We also 
plot the lightcone as a light-gray line as a 
visual aide. 
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FIG. 5: Parametric plot of the position of 
the unphysical /?-mode in the complex z plane 
for momentum k = {mD/'^,'mo,'^iTiD} show- 
ing that for k >^ rriu the unphysical /?-mode 
moves into the spacelike region of the real z-axis 
for large-^. 



Note, however, that the second term in Eq. (|2H) indicates the presence of unstable modes on the 
physical sheet for this angle of propagation and, as a result, the isotropic and anisotropic dispersion 
relations for the unphysical a-modes are not trivially connected. This is demonstrated by the fact 
that the type-1 mode (dotted lines) in Fig. do not extend down to A; = for finite ^ but instead 
terminate at a finite k = kg. Below ko the type-1 mode in the LHP moves onto the physical sheet 
and becomes the unstable (anti-damped) physical a- mode already discussed in Ref. [7|. The type-1 
mode in the UHP likewise moves onto the physical sheet and becomes the damped physical a-mode 
solution. 

A plot of the dispersion relations of the physical and unphysical /3-modes for ^ = 10 is shown 
in Fig. |1] (the mirror region Rez < is not shown). As one can see, the unphysical Plhp/uhp 
mode is lightlike for small momentum and spacelike at large momentum. As we will discuss below 
this mode is physically relevant if Rez is approximately spacelike, (Rez)^ < 1, and the modulus 
of Imz is small. In addition, in Fig. [S] we have plotted the position of the unphysical /3-mode pole 
for momentum k = {7n£)/2, mo, 2m£)} in the complex plane for various values of ^. 

From Fig. [H] we see that, as the anisotropy is decreased the Imz of the unphysical /3-mode 
becomes large at all momentum and therefore these unphysical modes will have a negligible impact 
on the propagator on the physical Riemann sheet. However, as the anisotropy is increased with 
fixed momentum, the modulus of Im z decreases so that these modes can become relevant for large 
anisotropies. From Fig. [S] we can also see that for fixed ^ and decreasing k that the unphysical 
/?- modes move to z = oo and so are unimportant based on the criteria stated above. However, 
for k > m£, the unphysical /3-modes have (Rez)^ ^ 1. Therefore for large anisotropies and k > 
rriD the effects of the unphysical /3-modes can have an important impact on the propagator for 
spacelike modes on the physical Riemann sheet. Additionally, we see that for infinite ^ the Imz of 
the unphysical /?-mode vanishes and it then moves to the physical sheet as we will discuss below. 
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FIG. 6: Sketch of the complex z plane in- 
cluding the extension of the logarithm to the 
unphysical sheet. Also shown is how a pole 
in the unphysical region (mountain) has ef- 
fects felt on the physical sheet. The black 
line indicates where the two sheets are joined 
together. 
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FIG. 7: Dispersion relations for 0„ = and 
^ = 10^: the /3-mode approaches its large-^ 
behavior. Again the physical a-, physical /3- 
, and unphysical /3-modes are indicated by 
dashed, solid, and dotted lines, respectively. 



C. Mountains on spirals 

As discussed in the previous section we find that, in addition to modes on the physical sheet 
found in Ref. [3|, for anisotropic systems there are also unphysical a- and /?-modes on neighboring 
Riemann sheets. But why bother about these modes, given that they do not "live" on the physical 
sheet? To see that there can be an effect on physical quantities, imagine the structure of the 
complex z plane spanned by the different sheets of the logarithm in the form of a spiral staircase: 
the physical sheet would correspond to the region covered by the spiral plane from the ground 
floor to the first floor, while the unphysical sheet where the extra quasiparticle mode lives would 
correspond to the region first to second floor. 

However, since for a range of momenta the unphysical /3-mode corresponds to an approximately 
spacelike pole, the propagator has a mountainous dent (singularity) in the spiral plane somewhere 
from the first to the second floor, with its peak the nearer the first floor the larger ^ is. But because 
the mountain has a finite width, its base can be felt also below the second floor, especially if its 
peak is near the first floor (see Fig. El for a sketch). This is precisely the situation we presented 
in previous section. There we showed that in the large-^ limit that for k ^ mo the Im z of the 
unphysical /3-modes becomes very small and (Re z)^ ;$ 1. Therefore, we expect that the unphysical 
modes on neighboring Riemann sheets do have physical consequences for large anisotropics, since 
then their effect can extend down to the physical sheet. In fact, as we will discuss below, in the 
limit of infinite ^ the pole moves onto the physical sheet itself. Note that the physical effect of any 
singularities existing on higher Riemann sheets (|n|> 1) would be negligible since the effect of these 
singularities would have to extend through the intermediate sheets prior to getting to the physical 
sheet, i.e. the base of the mountain would have to extend all the way down to the physical sheet 
through the spacelike spirals. 



D. Towards large ^ 

For very large values of ^, the dispersion relations are shown in Fig. [71 As one can see, the 
physical /3-mode now hits the lightcone at the point where the unphysical mode becomes spacelike. 
Indeed, when one takes the limit of ^ ^ oo, the /3-mode simply becomes 



1 
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0, 



(23) 



which has the simple real and propagating solution w^ = 7r?7i^/4. For any finite ^ the logarithmic 
singularity of (3 in Eq. (jlH) at the lightcone causes the physical /3-mode to always be timelike; 
however, the unphysical solution exists for both timelike and spacelike momenta. 



IV. SPECIAL CASE II : ^ 



oo 



Another special case where one can explicitly calculate the structure functions is when ^ — > oo. 
In this case it has been found in Ref. Isll that the distribution function becomes 



/oo . 

dx /iso {py/l + 
-oo ^ 



(24) 



which corresponds to the extreme anisotropic case considered by Arnold, Lenaghan, and Moore 
As a consequence, one can make use of this form by partially integrating Eq. (^ to obtain 



I?^{K) = 27ra, 



d^P Mp) 



(27r)3 p 



6'^ 



k^v^ + kH'^ {—up + k'^)v^v^ 



+ 



le 



(25) 



and applying the techniques from Ref. 
in the large-^ limit. Using 



lim 27rQs 

^^oo 



K-V-ie ' {-K-V 
to obtain analytic expressions for the structure functions 

d^p /c(p) 



(27r)3 p 



2 ^ 



(26) 



the structure functions are obtained using the contractions Eq. ©, giving 
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4 4 sin'' 
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6 + 2 cos Idn ± z 



3 - 6^2 - 2(1 + z^) cos 2Qn - cos 46l„ 
(z + sin6'„)3/2(2; - sin6'n)^/2 



,2\2 



z± 



(l-2z^)cOS^gn-(l 

(z + sin0„)3/2(^_sine„)3/2 



(27) 



where a positive sign above corresponds to the physical sheet and a negative sign corresponds to 
the unphysical sheet. Note that the lightcone singularity at z^ = 1 is not present in these structure 
functions any longer. In fact, all structure functions turn out to be purely real for z^ > sin^ 6 while 
there is a singularity located at z'^ = sin^ 6. Below this lightcone-like structure, the imaginary 
parts are non- vanishing. 
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FIG. 8: Dispersion relations for ^ = cxd for (a) 9n = j and (b) 0„ = §. Shown are the physical 
a— mode (short dashed lines) and the two physical /3-modes (full and dotted lines) as well as the 
unphysical /3-mode (long-dashed lines), respectively. For On = 0, the dispersion relations resemble 
those of Fig. [7] except that the mode which follows the light cone beyond fc^ = rriD'^/A in Fig. [7| 
ceases to exist and the /?-mode connects continuously across the lightcone to the now physical 
spacelike mode. 

A. Collective modes 

The dispersion relations of the collective modes are once more determined by the zeros of 
Eq. (|T!T|) . By conducting a Nyquist analysis similar to the previous section, one finds Na^phys = 
2 + 20 (— lim^^o ^a(-z)) and N^^^p^yg = 4 + 2Q {liuiz^o Ag{z)) for the physical sheet. These 
modes can be identified as the standard, propagating modes (one for Aq, and two for Aq together 
with their negative frequency equivalents) as well as two imaginary modes (one in the UHP and 
one in LHP) for Aq, and Aq, respectively. The latter again depend only on the sign of the static 
limit of the propagators, which in the limit of large ^ are functions of the momentum k and the 
angle On, respectively.^ 

By once again extending the structure functions Eq. H27|) to the unphysical sheets in the LHP 
and UHP (the terms involving square roots simply pick up an overall minus sign) we are also 
able to count the modes there by repeating the earlier analysis. One finds Na^LHP = Na,UHP = 
(lim^^o ^yl(^)) corresponding to a purely imaginary mode, and N/s.q,lhp = ^Aa,UHP = 2 — 
O (lim^^o ^g(^))- The first two of the latter modes are in general complex modes but may - 
in some restricted region of parameter space - also become purely imaginary solutions; the step 
function then just encodes the fact that one of the purely imaginary solutions moves to the physical 
sheet. 



V. CONCLUSIONS 

In this paper we have extended our studied the gluon polarization tensor in an anisotropic 
system. We extended our previous analysis to unphysical Riemann sheets and showed that for 
anisotropic distribution functions there are modes (singularities) in the "spacelike region" of the 
unphysical sheet which become physically relevant for large anisotropics. The chief way that these 



^ See Ref.jsl for a discussion of the static limit of Eq. if^ . 
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modes affect the physics is by altering the behavior of the propagator at soft spacehke momenta. 
The behavior of the propagator in this region determines the rate of energy transfer from soft 
to hard modes and the sign of this energy transfer may change for anisotropic systems so that 
there is instead a transfer of energy from hard to soft modes [igj. Whether or not the presence 
of the modes on the unphysical Riemann sheets are responsible for this "anti-Landau-damping" 



will be investigated in a separate paper |17[ in which we compute heavy fermion energy loss in a 
quark-gluon plasma along the lines of Ref. [l^ in which we calculated the same in an anisotropic 
QED plasma. 

In addition to these unphysical /3-modes we found that for anisotropic distributions the unphys- 
ical a-modes are different than in the isotropic case. In the isotropic case, there are two unphysical 
a-modes with one being in the upper half plane of the n = \ unphysical sheet and one in the 
lower half plane of the n = — 1 unphysical sheet. We showed that for finite-.^ and small momentum 
these isotropic unphysical modes move onto the physical sheet and become the unstable modes 
already discussed in Ref. ^7] . Additionally, for finite anisotropy there are two additional unphysical 
a-modes, again, with one being in the upper half plane and one in the lower half plane. These 
modes would, in principle, determine the dynamic screening of the magnetic gluon interaction at 
small momentum; however, the instabilities in this region will dominate the physics so it's not 
entirely obvious how important these additional unphysical a-modes are. 

In closing we would like to point out that although the finite-^ analysis presented here was 
performed only for On = there are generally relevant unphysical a- and /?-modes for all angles of 
propagation. The analysis proceeds exactly as discussed here but the details (number of modes etc.) 
changes here and there. Additionally, the results presented here are also applicable to anisotropic 
ultrarelativistic QED plasmas. 
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